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performed  in  ASWTDA  are  described,  followed  by  a  probabilistic  analysis.  In  the 
analysis  sections,  probabilistic  support  for  the  applied  methodology  is  provided  where 


applicable,  and  conceptual  problems  and  possible  solutions  are  cited  where 


iii 


TABLE  OF  CONTENTS 


I.  INTRODUCTION .  1 

A.  THESIS  GOAL .  1 

B.  ASWTDA  GOALS .  1 

C.  OVERVIEW  OF  ASWTDA  OPERATION  .  2 

D.  ASWTDA  BACKGROUND .  2 

E.  ORGANIZATION .  4 

II.  COMPUTER  ASSISTED  SEARCH,  A  CLASSICAL  MODEL .  5 

A.  PRIOR  PROBABILITY  MAP  .  5 

B.  TARGET  MOTION  MODEL  .  6 

C.  PROBABILITY  MAP  UPDATING  .  6 

1.  Search  Information  .  6 

2.  Target  Motion .  9 

III.  TERMINOLOGY  AND  NOTATION  .  10 

A.  TERMINOLOGY  .  10 

B.  NOTATION .  11 


iv 


IV.  DATA  FUSION  .  12 

A.  METHODOLOGY  .  12 

1.  Updating  Target  Density  for  Information .  12 

a.  Prior  to  Receipt  of  Positive  Information  (t<t0) .  12 

b.  Upon  Initial  Receipt  of  Positive  Information  (1=10) .  13 

2.  Subsequent  to  Receipt  of  Positive  Information  (t>to) .  15 

3.  Updating  Target  Density  for  Further  Positive  Information  ...  15 

B.  PROBABILISTIC  ANALYSIS  .  15 

V.  TARGET  MOTION  .  18 

A.  METHODOLOGY  .  18 

1.  Area  Clearance  Target  diffusion  .  18 

2.  Primary  Target  Diffusion .  21 

B.  PROBABILISTIC  ANALYSIS  .  23 

1.  Area  Clearance  Target  Diffusion .  23 

2.  Primary  Target  Diffusion .  31 

VI.  TARGET  COMBINATION  AND  CREDIBILITY .  33 

A.  TARGET  COMBINATION .  33 

1.  Single  Target .  33 

2.  Two  Separate  Targets .  35 

3.  Analysis  .  36 


a.  Single  Target  Model .  36 

b.  Two  Target  Model .  39 

4.  Modifications  to  Comply  With  a  One  or  Two  Target  Model  .  .  39 

B.  CREDIBILITY  .  41 

1.  A  Proper  Application  of  Credibility .  41 

2.  Credibility  in  ASWTDA .  42 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS  .  44 

LIST  OF  REFERENCES .  45 

INITIAL  DISTRIBUTION  LIST  .  46 


VI 


I.  INTRODUCTION 


A.  THESIS  GOAL 

The  goal  of  this  thesis  is  to  examine  the  methodology  used  in  ASWTDA  from 
a  probabilistic  standpoint.  It  is  intended  to  provide  probabilistic  support  for  the 
applied  methodology  where  applicable  and  to  cite  conceptual  problems  and  suggest 
possible  solutions  where  appropriate. 

B.  ASWTDA  GOALS 

ASWTDA  is  a  tool  to  assist  platform,  unit  or  force  commanders  afloat  and 
ashore  in  making  tactical  ASW  decisions.  It  is  necessary  that  it  support  rapid 
assimilation  of  data  pertinent  to  the  problem  and  needed  in  the  decision  making 
process.  ASWTDA  is  to  be  an  aid  in  the  Search  and  Detection  phases  of  the  ASW 
prosecution  effort  during  which  the  target’s  position  is  not  known  with  sufficient 
accuracy  for  weapon  placement  [Ref.  1:  p.  v].  Target  localization  and 
weapon  control  is  to  be  performed  by  existing,  onboard  fire  control  systems. 
ASWTDA  is  a  decision  support  tool  to  assist  in  the  conduct  of  the  following 
[Ref.  2:  p.  1]: 

1.  Large  area  environmental  analysis 

2.  Large  area  search  planning 
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3.  Resource  allocation  and  management 

4.  ASW  contact  correlation  and  management 


C.  OVERVIEW  OF  ASWTDA  OPERATION 

ASWTDA  begins  operation  at  some  problem  start  time  with  an  initial  target 
location  distribution  and  an  assumed  target  motion  model.  ASWTDA  then  processes 
data  from  negative  search  efforts,  positive  contact  reports  and  target  motion  to 
produce,  at  some  future  forecast  time,  a  probability  distribution  for  the  target’s 
position.  This  data  may  then  be  displayed  as  a  probability  map  showing  regions 
where  target  location  is  more  or  less  likely,  or  it  may  be  combined  with 
environmental  and  sensor  data  to  produce  a  map  showing  regions  where  target 
detection  is  more  or  less  likely. 

D.  ASWTDA  BACKGROUND 

The  development  of  ASWTDA  began  about  November  1988.  Chief  of  Naval 
Operations  OP-71,  Antisubmarine  Warfare  Division,  was  named  sponsor  of  the 
program,  NAVSEA  63-D3  was  named  program  manager.  OP-71  and  NAVSEA  then 
named  Naval  Oceanographic  Systems  Center  (NOSC)  of  San  Diego,  CA  as  Technical 
Director,  Naval  Underwater  Systems  Center  (NUSC),  New  London,  CT  as  developer 
and  Naval  Air  Development  Center  (NA DC),  Warminster,  PA  as  the  laboratory  for 
independent  verification  and  validation.  [Ref.  3] 
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It  was  felt  that  existing  computer  search  aids  were  not  adequate  for  the  ASW 
commanders  and  operators  at  sea.  OP-71  established  a  fleet  working  group  to 
provide  requirements  for  a  new  ASWTDA  to  meet  the  Navy’s  needs.  A  major  goal 
was  to  get  the  program  to  the  fleet  as  quickly  as  possible.  The  development  was  to 
be  an  evolving  process  with  rapid  prototyping  and  fleet  feedback.  In  order  to  reduce 
development  effort,  the  program  was  to  evolve  from  three  previous  efforts.  These 
were  [Ref.  3]: 

1.  Integrated  Tactical  Decision  Aid  (ITDA) 

2.  Composite  Area  Search  Evaluation  (CASE) 

3.  Fast  Frigate  Integrated  Sonar  Tactical  System  (FFISTS). 

The  analysis  performed  in  this  thesis  is  based  on  the  methodology  as  described 
in  Reference  1.  This  was  the  methodology  in  use  as  of  November  1989.  Due  to  the 
nature  of  the  development  process  of  ASWTDA,  the  program  may  change  greatly  in 
a  short  period  of  time.  The  reader  must  bear  in  mind  that  this  is  based  upon  dated 
information  and  the  program  as  analyzed  is  not  necessarily  what  is  being,  or  has  been 
delivered  to  the  Navy  for  fleet  use.  As  of  this  writing,  a  version  of  ASWTDA  is  in 
place  for  at  sea  testing  onboard  the  U.S.S.  Cushing,  the  flagship  of  Destroyer 
Squadron  31,  an  ASW  squadron  based  in  San  Diego,  California. 


E.  ORGANIZATION 


Chapter  II  describes  a  model  of  a  classical  Computer  Aided  Search  (CAS) 
program.  This  will  provide  the  reader  with  a  general  background  of  the  operations 
performed  in  such  a  program  and  the  notation  used  in  the  description  of  the 
operations  of  ASWTDA. 

In  Chapters  IV  and  V,  the  following  specific  areas  will  be  addressed  in  this 
analysis: 

1.  Data  Fusion 

2.  Target  Motion 

Each  chapter  is  broken  down  into  two  major  sections.  In  the  first  section, 
Methodology,  a  description  is  given  of  the  mechanics  of  the  operations  involved,  with 
no  attempt  to  discuss  the  reasoning  behind  the  operations.  In  the  second  section, 
Probabilistic  Analysis  an  attempt  is  made  to  provide  probabilistic  justification  for  the 
methods  employed. 
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II.  COMPUTER  ASSISTED  SEARCH,  A  CLASSICAL  MODEL 


A  Classical  model  of  a  moving  target  CAS  program  is  a  program  which  outputs 
a  probability  map  of  target  position  at  any  time  specified  by  the  operator.  The  three 
main  elements  of  a  Markoviar  moving  target  CAS  system  are: 

1.  Prior  Probability  map  of  target  Position 

2.  Model  of  Target  Motion 

3.  Updates  of  the  probability  map  for  target  motion  and  information  received  on 
the  target 

Advanced  programs  of  this  type  may  also  provide  the  user  with  search  plan 
recommendations. 

A.  PRIOR  PROBABILITY  MAP 

In  CAS  applications,  the  geographic  region  of  interest  is  always  divided  into  an 
array  of  cells.  The  target  location  is  described  by  a  probability  distribution  over  these 
cells.  We  will  index  the  cells  in  the  array  with  the  subscripts  i,j:  i=l...N,  j=  1...M;  this 
set  of  cells  makes  up  the  analysis  aiea,  5.  Let  E;j(t)  be  the  event  a  target  is  in  cell 
(i,j)  at  time  (t)  and  define  pjj(0)  =  P(Eij(0)),  ie.  p^O)  is  the  probability  the  target  is 
located  in  cell  (i,j)  at  time  (0). 


The  prior  probability  map  is  constructed  by  assigning  a  number  between  0  and 
1  to  each  ptj  at  time  t=0  such  that: 

N  M 

EEtf)'10  (1) 

i*l  j* 1 

This  distribution  is  called  the  "prior  probability  map  of  the  target",  frequently  referred 
to  simply  as  the  "prior".  If  the  operator  has  no  information  about  target  location,  the 
prior  may  be  simply  a  uniform  target  distribution. 

B.  TARGET  MOTION  MODEL 

Target  motion  must  be  described  in  probabilistic  terms  because  of  the 
uncertainties  involved  with  the  motion.  A  common  motion  model  that  is  used  is  the 
Markov-chain  model.  The  target  location  is  known  as  the  state  of  the  target.  Some 
proportion  of  the  probability  moves  from  each  cell  to  the  others  in  each  time  step, 
as  determined  by  a  transition  matrix,  Q.  The  Markov-chain  model  contains  no  target 
history  other  than  what  is  contained  in  the  state.  Thus  the  cell  probabilities  change 
at  a  given  time  based  only  on  the  information  available  in  the  cell  at  that  time  (the 
definition  of  a  Markov  process)  [Ref.  4:  p.  21]. 

C.  PROBABILITY  MAP  UPDATING 
1.  Search  Information 

Search  information  may  be  of  two  types:  1.  Negative  information:  search 
is  performed,  but  no  target  was  found.  2.  Positive  information:  a  report  on  target 
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position  resulted  from  search.  For  our  purposes,  local  sensors  (sonobuoys,  towed 
arrays,  etc.)  provide  only  negative  information.  If  target  detection  occurs  by  these 
sensors,  the  problem  is  one  of  localization  and  is  no  longer  a  wide  area  search 
problem  of  the  type  addressed  by  ASWTDA.  Positive  information  may  be  obtained 
from  a  long  range  type  sensor,  e.g.  SOSUS,  which  does  not  provide  a  target  position 
with  sufficient  accuracy  for  target  localization  without  further  search. 

Updating  for  search  information  is  done  by  application  of  Bayes’  theorem. 
Let  I(t)  be  the  event  that  some  information  has  been  received  about  the  target  at 
time  (t).  Assuming  that  P(I(t)  |  Ejj(t)),  i.e.,  the  probability  that  information  I(t)  is 
received,  given  the  target  is  in  cell  (i,j)  at  time  (t),  is  known  for  all  i,j,  the  target’s 
position  distribution  can  be  modified  on  the  basis  of  receipt  of  this  information. 
These  probabilities  are  based  on  models  of  search  effectiveness  for  all  types  of  search 
performed  in  the  analysis  area.  The  target  distribution  is  modified  as  follows 
[Ref.  5:  p.  65]: 


p('(»  |g|,<»)W>) 
E  '’(«<)  I5j«)W>) 

iJeS 


(2) 


Note  that  P(Ej  j(t)  |  I(t))=P(Ej  j(t))  when  P(I(t)  |  E4 j(t))  is  independent  of  (i,j).  This 
is  commonly  assumed  to  be  the  case  when  I(t)  is  "null",  i.e.,  no  information  is 
received. 


For  the  purposes  of  this  paper,  information  is  received  about  the  target  at 
discrete  times  At  apart.  These  times  will  be  denoted  by  t,t+l,t+2,....  Let: 
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Equation  (2)  becomes: 


pjt)  =  P(Ei/t)\I(l)f.J(t-l)) 

Pi/t)  = 


(3) 


P(/(Q  1  EJt))PiJt) 

E  *>('(')  IW> 

iJtS 


(4) 


This  requires  an  assumption  about  independence,  however.  Note  that  the  p,j(t)  are 
actually  conditional  probabilities,  with  the  condition  being  the  information  history 
H,=(I(0),...,I(t-l)),  prior  to  receipt  of  I(t).  Thus  pjj(t)=P(Eij(t)  |  Ht),  and  the  problem 
is  to  calculate  P(Ejj(t)  |  I(t)H().  Bayes’  theorem  is  actually: 


P(Etj(t)  |  H'I(t)) 


P(I(t)  1  Ht  EJtmE^t)  1  Ht)P(Ht) 
52  P(I(t) \HtEi/t))P(Eij(t)  I H')P{Ht) 


(5) 


where  P(HlE1J(t))=P(EiJ(t)  |  H,)P(H,).  We  can  then  cancel  P(H,)  and  Equation  (5) 
becomes: 


P(EJt)\H,I(t)) 


P(7(Q  (0)^(0 

j:pm\HtEjmpjf> 


(6) 


which  is  equivalent  to  Equation  (2)  if  it  is  true  that  P(I(t)  { HtEjj(t))= P(I(t)  ( E;j(t)); 
i.e.,  conditional  on  target  location  being  given,  I(t)  is  independent  of  H,.  This 


assumption  is  implicit  in  the  repeated  application  of  Bayes’  theorem,  i.e.,  each  new 


piece  of  information  must  be  conditionally  independent  of  past  information.  This  is 
usually  a  reasonable  assumption,  with  the  exception  of  systems  such  as  SOSUS,  which 
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gives  frequent  reports  that  are  correlated,  even  when  target  position  is  given. 
[Ref.  6:  p.  7] 

2.  Target  Motion 

In  the  Markov-chain  model,  the  motion  update  requires  an  update  of  the 
target  probabilities  in  each  cell.  Given  that  a  target  is  in  cell  (i,j),  it  may,  in  a  given 
time  interval,  remain  in  cell  (i,j),  or  transition  to  another  cell  in  the  analysis  area. 
Associated  with  each  cell  (i,j),  is  a  transition  matrix,  Q,  which  describes  target  motion 
from  cell  (i,j)  in  one  time  step.  Let  Q(i,jjm,n)  be  P(Ejj(t+l)  |  Emn(t)).  Then, 
assuming  Markov  target  motion  [Ref.  4:  p.  240]: 

PiP*  1)  =  E  E  Pm.nW  07  I  »»,«)  (7) 

neS  m  eS 

Thus,  as  the  notation  is  used  here,  p^t+l)  is  obtained  from  pT^t)  by 
application  of  the  transition  matrix,  Q,  for  target  motion  during  the  time  interval  (t) 
to  (t+1).  P+,j(0  is  obtained  from  p,j(t)  by  application  of  Bayes’  theorem  for 

information  received  at  time  (t).  The  sequence  of  computations  thus  goes:  p(0), 
p+(0),  p(l),  p+(l),  p(2),  etc. 

In  principle,  the  target  state  maybe  much  more  complex  than  simply  target 
location.  It  may  contain  other  properties  such  as  target  course,  depth,  motion  model, 
time,  etc.  The  transition  matrix,  Q,  would  then  contain  transition  probabilities  for 
each  of  these  state:;  to  every  other  possible  state.  This  more  general  notation  will  not 
be  used  in  this  thesis  however,  because  the  state  in  ASWTDA  consists  only  of 
position. 
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III.  TERMINOLOGY  AND  NOTATION 

A.  TERMINOLOGY 

ASWTDA  performs  analysis  within  a  geographic  region  with  location  and  size 
specified  by  the  operator  at  the  beginning  of  the  session.  The  region  is  a  square 
ranging  from  100  to  500  Nautical  Miles  on  a  side.  This  region,  known  as  the  Analysis 
Area,  is  then  fixed  and  may  not  be  changed  without  problem  re-initialization.  This 
analysis  area  is  subdivided  into  a  grid  of  1764  (42  by  42)  square  cells.  Thus,  cell  size 
dimensions  range  from  2.38  to  11.90  nautical  miles.  Target  location  is  maintained  in 
ASWTDA  as  a  probability  distribution  over  the  cells. 

ASWTDA  maintains  two  target  density  matrices: 

1.  Area  Clearance  Density  Matrix 

2.  Primary  Target  Density  Matrix 

The  purpose  of  the  Area  Clearance  target  density  display  is  to  illustrate  target 
density  in  the  absence  of  target  detection.  It  is  a  display  to  be  used  primarily  for 
determining  regions  which  have  been  searched  and  sanitized,  thereby  containing  a  low 
probability  of  a  submarine  threat.  These  are  areas  into  which  the  operational 
commander  may  safely  send  vulnerable  forces. 

The  Primary  Target  density  display  shows  the  Primary  Target  density  at  any 
time  after  target  detection  has  occurred.  This  is  intended  to  be  the  display  to  which 
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the  operator  focuses  his  attention  following  target  detection,  although  the  Area 
Clearance  display  may  be  monitored  for  the  reasons  stated  above. 

The  operator  may  also  view  a  combined  density  display  which  is  formed  from 
the  sum  of  the  Area  Clearance  and  Primary  Target  density  matrices  and  is  calculated 
at  the  time  of  display. 

B.  NOTATION 

Each  density  matrix  is  described  as  a  42  by  42  cell  array  with  AQj  representing 
the  probability  density  for  a  particular  cell  in  the  Area  Clearance  matrix,  and  PTSj  the 
probability  density  for  a  cell  of  the  Primary  Target  density  matrix. 

For  a  matrix  X,  the  norm,  |  X  | ,  is  defined  by: 

*1  =  £  X.J  (8) 

Also,  the  product  XxY  for  matrices  X  and  Y  is  a  multiplication  on  a  cell  by  cell  basis 
of  the  values  of  X,  j  and  Yij7  rather  than  the  matrix  product  of  Linear  Algebra. 

The  notation  as  described  in  Chapter  II  with  regard  to  information  and  position 
updating  will  be  utilized  for  describing  the  operations  through  time  on  the  target 
density  matrices. 

For  any  two  events  x  and  y,  when  these  are  listed  together  in  a  probability 
statement,  the  meaning  is  (x  and  y),  e.g.,  P(xy)=P(x  and  y). 
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IV.  DATA  FUSION 


A.  METHODOLOGY 

The  target  density  matrices  are  updated  for  information  at  6-minute  intervals. 
This  update  occurs  at  the  end  of  each  6-minute  period.  At  each  update,  the  negative 
effects  of  all  search  efforts  are  incorporated,  then  diffusion  due  to  target  motion  is 
applied. 

ASWTDA  incorporates  negative  information  from  three  types  of  acoustic 
search:  sonobuoys  (single  sonobuoys  and  sonobuoy  fields),  search  tracks  and  search 
areas.  A  search  track  is  used  when  the  sensor  platform’s  track  is  known  with 
reasonable  accuracy.  A  search  area  is  used  when  a  bounded  region  is  assigned  to  a 
search  platform,  but  the  platform’s  motion  within  the  region  is  not  specified. 
[Ref.  1:  p.  53] 

1.  Updating  Target  Density  for  Information 

a.  Prior  to  Receipt  of  Positive  Information  (t<to) 

On  problem  initialization,  ASWTDA  assumes  a  single  target  is 
distributed  uniformly  over  the  analysis  area,  thus  Pij(0)=  1/N  for  all  non-land  cells, 
where  N  is  the  number  of  non-land  cells  within  the  analysis  area,  S.  Prior  to  receipt 
of  positive  information  on  a  target,  this  probability  distribution  is  contained  entirely 
within  the  Area  Clearance  target  density  matrix,  thus  ACjj(0)=Pjj(0)  for  all  (i,j)eS, 
and  the  Primary  Target  density  matrix  is  initialized  to  0,  so,  PTjj(0)=0  for  all  (i,j)eS. 
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Let  I(t)  be  the  event  that  no  target  detection  occurs  between  times 
(t-1)  and  (t).  For  a  position  update  at  time  (t),  the  cumulative  probability  of  detection 
for  each  cell  over  the  six-minute  time  interval  from  (t-1)  to  (t)  is  calculated  based  on 
all  search  efforts  in  progress  during  the  time  interval  from  (t-1)  to  (t).  These 
calculations  are  based  on  search  models  for  each  type  of  search  and  incorporate 
environmental  data  as  well  as  data  about  the  sensors  themselves.  Let 
Pnjj(t)=P(I(t)  |  Ejj(t)),  i.e.,  Pn^t)  is  the  probability  that  no  detection  occurs  in  cell 
(i,j)  during  the  time  interval  (t-1)  to  (t),  given  the  target  is  in  cell  (i,j).  The  Area 
Clearance  target  density  matrix  is  then  updated  by  applying  Bayes  theorem  to  all  cells 
as  follows: 

AC'(t)  =  ACM**#)  (») 

\AC(t)<P«(t)\ 

b.  Upon  Initial  Receipt  of  Positive  Information  (t—tff 

ASWTDA  assumes  that  the  target  distribution  based  on  the  positive 
information  can  be  described  by  a  bivariate  normal  distribution.  Let  B  represent  this 
bivariate  normal  distribution.  This  distribution  is  discretized  into  the  same  size  cells 
as  used  in  the  ASWTDA  analysis  area  and  probability  values  are  assigned  to  the  cells 
according  to  the  amount  of  the  bivariate  normal  density  encompassed  by  the  cell. 
The  discretized  bivariate  normal  density  will  be  denoted  B°.  Note  that  B°  extends 
infinitely  in  all  directions,  whereas  5  covers  only  a  finite  area. 

This  target  distribution  is  assigned  a  confidence  factor  (c/),  currently 
set  at  a  value  of  0.8  but  designed  to  be  an  operator  assigned  vaiue  in  future  program 
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releases.  The  discretized  bivariate  normal  distribution  is  renormalized  so  that 
|  B°  j  =cf.  The  confidence  factor  is  intended  to  be  representative  of  the  input 
distribution  as  a  whole  and  not  a  function  of  any  individual  contact  reports  used  to 
generate  the  input  distribution.  [Ref.  1:  p.  74] 

As  a  result  of  this  positive  information,  PTij(t0)=fiZ3ij  for  all  (i,j)eS, 
provided  the  cell  does  not  fall  on  land,  in  which  case  it  is  assigned  a  value  of  zero. 
The  distribution  contained  in 
PT(t0)  then  is  that  subset  of  B° 
which  is  contained  in  S  as 
illustrated  in  Figure  1.  Let 
pr=|PT(t0)|.  pr  will  be 
dependent  upon  the  amount  of 
Bd  which  falls  within  5,  not  on  a 
land  mass.  In  all  cases,  pr<cf ; 
for  example,  for  the  case 
illustrated  in  Figure  l,pnzcf/ 2. 

The  following 
operation  is  then  performed  to 

fuse  the  previous  negative  information  with  this  positive  information: 


Figure  1  Bivariate  Normal  Distribution  Overlaid 
on  Analysis  Area,  S. 


PT*(tj  = 


prxI^C+itJxPnt,))) 

|^C^(f0)xP71[r0))| 


(10) 
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2.  Subsequent  to  Receipt  of  Positive  Information  (t>t0) 

For  each  six-minute  time  interval,  negative  search  information  is 
incorporated  into  the  Area  Clearance  density  matrix  using  Equation  (9)  and  the 
Primary  Target  density  matrix  using  Bayesian  updating  in  the  same  manner  as 
follows: 


pr(t)  = 


PT{t)xPn(f) 
\P7\t)xPn(t)  | 


(ID 


3.  Updating  Target  Density  for  Further  Positive  Information 

Currently,  additional  positive  information  may  be  treated  in  one  of  two 


ways: 


1.  The  new  report  may  be  incorporated  into  the  Over  the  Horizon-Maneuvering 
Target  Statistical  Tracker  (OTH-MTST)  input  data  set  to  produce  a  new 
initial  target  distribution. 

2.  The  contact  report  may  be  used  as  the  goal  of  a  constrained  random  walk 
scenario,  with  a  past  report  as  the  initial  distribution  for  the  scenario.  In  this 
manner,  several  reports  may  be  incorporated  in  a  series  of  constrained 
random  walks. 


B.  PROBABILISTIC  ANALYSIS 

As  described  above,  negative  information  is  incorporated  into  both  target 
density  matrices  according  to  Bayes’  theorem.  It  can  be  shown  that  the  initia1  positive 
information  is  also  incorporated  in  a  Bayesian  manner  if  certain  assumptions  are 
made  about  the  received  information,  and  the  target  is  assumed  to  have  followed  the 
motion  model  of  the  Area  Clearance  target  prior  to  the  information  update. 
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Specifically,  if  at  time  (t0),  it  is  true  that  the  distribution  of  the  Primary  Target  is 
equivalent  to  the  distribution  of  the  Area  Clearance  target,  and  pr=  1,  then  Equation 
(10)  is  a  proper  application  of  Bayes’  theorem. 

There  are  two  issues  that  arise,  however: 

1.  Why  is  the  distribution  of  the  Primary  target  equivalent  to  the  distribution  of 
the  Area  Clearance  target  prior  to  the  positive  information  update? 

2.  Why  is  pr*  1? 

These  issues  will  be  discussed  in  Chapter  VI. 

Further  positive  information,  however,  is  not  incorporated  in  a  Bayesian 
manner.  Either  of  the  methods  mentioned  above  which  may  be  utilized  to 
incorporate  this  further  positive  information  results  in  a  change  of  the  Primary  Target 
motion  model. 

It  is  possible  that,  given  more  than  one  possible  target  motion  model,  receipt 
of  positive  information  may  alter  the  probabilities  that  we  assign  to  each  of  these 
models.  In  other  words,  we  change  our  belief  as  to  what  the  Primary  Target  is,  and 
has  been  doing,  based  on  this  information.  This  is  an  important  point  in  that  the 
target  motion  prior  to  the  receipt  of  the  positive  information  must  also  be  taken  into 
account.  This  is  not  the  case  in  ASWTDA.  Based  on  receipt  of  positive  information 
at  time  (t0),  the  Primary  Target  changes  its  motion  model,  i.e.,  at  time  (to),  the 
Primary  Target  distribution  is  dependent  upon  the  Area  Clearance  target  motion 
model  prior  to  time  (t0).  For  t>t0,  the  Primary  Target  then  follows  an  entirely 
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different  motion  model.  The  Primary  Target  can  not  be  expected  to  actually  alter 
its  motion  based  solely  on  the  receipt  of  information  on  the  target,  unless  it  is 
believed  that  the  target  is  reacting  to  our  search  efforts. 
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V.  TARGET  MOTION 


A.  METHODOLOGY 

ASWTDA  models  target  motion  using  a  Markov-Chain  approach.  Note  that 
in  the  general  notation  of  Markov  motion  as  described  in  Section  II.C.2,  the  transition 
matrix,  Q,  is  a  1764  by  1764  matrix,  i.e.,  it  describes  the  probability  of  transition  from 
any  state  to  any  other  possible  state.  In  ASWTDA,  however,  the  target  is  limited  to 
traveling  at  most  one  cell  length.  This  requires  a  much  smaller  transition  matrix, 
specifically  a  3x3  transition  matrix,  q^,  for  each  cell  (i,j)eS.  This  transition  matrix 
describes  target  motion  from  that  cell  during  any  single  motion  time  interval.  A 
portion  of  the  cell’s  target  density  remains  in  the  original  cell;  the  rest  is  distributed 
to  the  adjacent  eight  cells.  [Ref.  1:  p.  23] 

1.  Area  Clearance  Target  diffusion 

ASWTDA  assumes  that  a  target  associated  with  the  Area  Clearance 
Display  moves  about  randomly.  The  transition  matrix,  qAC,  associated  with  this 
motion  is  shown  in  Figure  2.  This  transition  matrix  is  applied  to  every  cell  within  the 
analysis  area,  thus  the  lack  of  the  i,j  subscript.  The  transition  matrix  coefficients  are 
fixed,  regardless  of  cell  size  or  target  speed. 

The  transition  matrix  coefficients  were  calculated  as  follows:  [Ref.  1:  p.  26] 

1.  Build  a  3x3  matrix  of  9  square  cells. 
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2.  Populate  the  center  cell  with 
1,000,000  uniformly  distributed 
points. 

3.  Move  each  point  one  cell-length 
in  a  random  direction  (from  0  to 
360  degrees). 

4.  Count  the  number  of  points  in 
each  of  the  9  cells. 

5.  Normalize  the  final  distribution 
of  points  and  adjust  to  ensure 
symmetry. 

At  each  diffusion  iteration, 
the  target  density  is  re-computed  using 
a  temporary  density  matrix.  Let  the 
temporary  matrix  be  called  T.  The  transition  matrix  is  applied  as  follows: 


b 

c 

b 

c 

a 

c 

b 

c 

b 

a=4.4939999E-2 
b=7.9482503E-2 
c-1 .5928250E-1 

Figure  2  Area  clearance  Transition  Matrix 


1.  Zero  the  temporary  matrix,  1. 

2.  For  a  cell  (i,j)  of  the  Area  Clearance  density  matrix,  AC,  multiply  each  of  the 
values  of  the  transition  matrix,  q,  by  AC^. 

3.  Add  these  values  to  the  matrix  T.  Thus,  the  density  from  cell  (i,j)  is 
distributed  in  cell  (i,j)  and  its  eight  surrounding  cells  in  matrix  T.  (See 
Figure  3) 

4.  Repeat  for  every  cell  (i.j). 

5.  Replace  AC  with  T. 

6.  Repeat  the  process  if  more  transition  matrix  applications  are  required. 
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At  edge  and  corner  cells,  any  density 
that  would  diffuse  out  of  the  analysis  area  is 
retained  in  the  original  cell  (i,j).  At  land 
boundaries,  any  density  that  would  diffuse 
onto  land  is  distributed  into  cell  (i,j)  and  its 
immediate  surrounding  non-land  cells.  Note 
that  analysis  area  boundary  and  land  effects 
do  not  change  the  transition  matrix  qAC,  but 
are  accounted  for  at  each  application  of  qAC. 

We  will  use  the  following  compact  notation 
to  represent  this: 

p(t*  1)  =  p*(t)+q  (*2) 

Because  the  coefficients  for  this  matrix  are  fixed,  the  Area  Clearance  target 
transition  matrix  is  applied  using  Equation  (12),  an  irregular  number  of  times  per  six- 
minute  time  interval,  determined  by  target  speed,  mean  time  on  a  given  course  and 
cell  size.  In  this  discussion,  an  iteration  is  one  six-minute  interval,  while  an  application 
of  the  transition  matrix  refers  to  a  single  use  of  equation  (12).  This  update  occurs 
as  follows  [Ref.  1:  p.  26]: 

1.  Calculate  R86,  the  radius,  in  nautical  miles,  which  contains  86  percent  of  the 
target  density.  For  a  given  target  speed  (V)  in  Kts.,  mean  time  on  a  given 
course  (X)  in  hours  and  length  of  time  (T)  in  hours,  R86  is  approximated  by: 


— 

AC,, 
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cACj  | 
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cACi.i 

bAQ 

a* 

4b*4 

C*1 

Figure  3  Application  of  Area 
Clearance  Transition  Matrix 
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R86  =  2„ 


,  v*\2 

f  _r  A 

e  A-l  +  - 

N 

V  ^  / 

min(1.0 , 0.707  +  0.293  T/  3  A) 


(13) 


2.  Calculate  the  number  of  transition  matrix  applications,  n,  that  would  produce 
an  R86  equal  to  the  R86  approximated  by  the  calculation  in  step  1,  where: 

=  (R86/Ce  11  Length)2  (14) 

2.45 

3.  Since  n  may  not  be  an  integer  value,  the  calculated  value  is  rounded  to  the 
nearest  whole  number.  If  the  required  number  of  transition  matrix 
applications  for  the  current  problem  time  (total  time  from  t=0  to  time  t  )is 
greater  than  the  number  of  applications  performed  as  of  the  end  of  the 
previous  iteration,  then  the  diffusion  process  is  called  and  the  transition  matrix 
is  repeatedly  applied  using  Equation  (12)  until  the  number  of  applications 
performed  equals  the  required  number  of  applications  computed  in  step  2. 


2.  Primary  Target  Diffusion 

For  Primary  Target  displays,  a  unique  3x3  transition  matrix  is  calculated  for 
each  cell  within  the  analysis  area.  These  transition  matrices  will  be  denoted  by  qPTjj. 
These  matrices  are  calculated  by  performing  a  Monte  Carlo  simulation  using  2000 
tracks  originating  from  a  specified  initial  target  distribution,  and  proceeding  through 
an  anticipated  target  motion  scenario.  Land  avoidance  is  performed  in  the  process 
of  the  Monte  Carlo  simulation  thus  preventing  the  Primary  Target  from  diffusing 
onto  land.  As  a  given  Monte  Carlo  track  proceeds,  the  transition  matrix  building 
subroutine  keeps  track  of  the  number  of  transitions  from  any  one  cell  to  an  adjacent 
cell  in  the  matrix  at  intervals  of  dt,  given  by: 
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(15) 


(Average  Target  SpeedXCos(6)) 


where  <5  is  the  difference  from 
the  average  target  course  to  the 
nearest  cardinal  heading,  so 
-7r/4  <  <5  <  tt/4.  The  cosine  term  is 
used  because  it  was  determined 
experimentally  that  a  target  on  a 
diagonal  course  tended  to  diffuse 
too  quickly  in  actual  application. 

The  simulation  process  is 
illustrated  in  Figure  4.  There  is 
a  3x3  matrix  associated  with  each 
of  the  cells  in  the  upper  matrix. 

The  target  track  is  marked  off  in 
numbered  intervals  of  dt.  Thus 
after  completion  of  this  track,  the  3x3  matrix  for  cell  (1,1)  would  contain  the  numbers 
shown  in  the  matrix  at  the  bottom.  Because  dt  is  calculated  using  the  average  target 
speed,  it  is  possible  that  the  target  may  travel  farther  than  one  cell  length  in  a  time 
step  dt,  this  is  illustrated  in  time  step  8  of  Figure  4.  If  this  occurs,  a  count  is  placed 
in  the  nearest  adjacent  cell.  The  simulation  then  continues  from  the  actual  simulated 
position.  After  all  2000  tracks  have  been  generated,  the  cell  counts  for  each  3x3 


Figure  4  Simulation  of  One  Track  for  Primary 
Target  Diffusion 
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matrix  are  translated  into  transition  probabilities.  This  is  performed  by  replacing 
each  cell  of  the  3x3  transition  matrix  with  the  number  of  counts  in  the  cell  divided 
by  the  total  number  of  counts  in  the  3x3  matrix. 

This  set  of  transition  matrices  is  used  until  additional  positive  information 
is  received  or  the  target  motion  scenario  is  modified,  at  which  time  a  new  Monte- 
Carlo  simulation  is  performed  and  the  transition  matrices  are  re-computed.  The 
Primary  Target  transition  matrices,  qPTy,  are  applied  using  Equation  (12)  at  regular 
time  intervals  equal  to  dt  as  calculated  in  Equation  (15). 

B.  PROBABILISTIC  ANALYSIS 

1.  Area  Clearance  Target  Diffusion 

The  model  that  ASWTDA  assumes  for  the  Area  Clearance  target  is  that 
of  a  Random  Tour.  Under  this  model,  the  target  is  assumed  to  move  at  constant 
speed,  but  to  make  course  changes  at  random  times.  Each  new  course  is  chosen 
from  a  uniform  distribution  on  [0,2tt].  The  random  travel  time  between  course 
changes  is  assumed  to  be  drawn  from  an  exponential  distribution.  Each  new  course 
and  time  on  course  are  assumed  independent.  Note  that  this  is  not  the  same  as  a 
diffusion  process,  which  is  not  a  constant  velocity  process. 

ASWTDA  uses  the  results  for  a  diffusion  process  to  approximate  the 
random  tour.  This  approximation  is  based  on  the  fact  that  the  variance  of  the 
diffusion  process  has  the  same  functional  form  as  E[R2T],  the  mean  squared  radial 
displacement  of  a  randomly  touring  target  [Ref.  7]: 
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(16) 


E  [/#  =  2V2X2 


Furthermore,  for  the  diffusion  model,  the  target  at  all  times  has  a  circular  normal 
distribution  with  mean  /j,=0  and  a2=1/2E[R2T].  For  a  circular  normal  distribution, 
R2j/2a2  has  a  unit  exponential  distribution  for  which  86  percent  of  the  density  is 
contained  within  2a.  [Ref.  5:  p.  312].  For  the  random  tour  model,  the  normal 
assumption  holds  only  for  large  T.  Therefore  the  correction  factor  for  T/A<3.0  is 
applied  in  equation  (13). 

It  is  possible,  however,  to  derive  an  exact  expression  for  R86.  This  is  done 
by  integrating  the  probability  density  function  of  the  particle’s  position  at  time  t,  given 
by  the  formula: 


,  -{(i-v'TV) 

g(W)  =  - *  for  p<l  (17) 

2tiA.(V7)V1-p2 

where  ps(x2+y2)/(VT)2  [Ref.  8],  Setting  up  the  integral  and  converting  to 
polar  coordinates,  for  z<VT  we  have: 
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P(Rr<.z ) 


/ 


u  T 

yjl-U1  * 


-I(i-vTv) 
x  du 


Substituting  the  following: 


v  =  y/l-U1 
dv  =  — ---  du 

-U1 


(19) 


(20) 


we  get: 


P(Rt<lz)  = 


=  e 


T 

A. 


l-e  x 


(21) 


=  l 

The  discontinuity  at  ^  =  1  is  due  to  the  fact  that  with  probability  e'm  the  target  does 
not  change  course,  so  there  is  a  probability  of  mass  e‘TO  concentrated  on  the  circle 
of  radius  VT. 

We  now  solve  for  R862rt,  the  radius  which  contains  86  percent  of  the 
target  density  for  the  random  tour  case: 
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(22) 


where  the  integer  2  appears  because  -ln(  1-0.86) =2.  Therefore: 

R86lT  =  V2!2  1-fl-— T 

"  l  t)\  (24) 

=  4V2  XT - 4V2  X2  \Tz2X 

R862rt  =  V*T2  ,T<2X 

In  Order  to  relate  R86  and  the  number  of  diffusion  iterations  performed 
for  the  Area  Clearance  target,  the  developers  of  ASWTDA  performed  a  series  of 
experiments  by  applying  the  Area  Clearance  transition  matrix  repeatedly  and  fitting 
the  data  experimentally,  obtaining  the  relation  of  Equation  (14).  Equation  (14)  can 
be  derived  analytically  however.  Furthermore,  it  can  be  shown  that  the  values  for  a,b 
and  c  in  qAC  can  be  varied  so  that  the  number  of  iterations  required  provides  for 
more  convenient  application  of  the  transition  matrix. 

First,  we  must  show  that  the  application  of  Markov  target  motion  used  in 
ASWTDA  produces  an  approximately  circular  distribution.  Assume  we  have  a  target 
at  position  (0,0).  Let  X„  be  the  horizontal  position  at  time  n.  Then: 
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(25) 


x.  -  x>. 

1*1 

where  u;  is  the  horizontal  movement  at  time  i.  In  ASWTDA,  since  the  target  may 
move  at  most  one  cell  length,  UjC(-l,0,l)  and,  due  to  the  symmetry  of  qAC  its  expected 
value  E[Uj]=0.  The  variance,  Var[uJ  =  E[u2i]-E[ui]2=E[u2i].  Using  the  values  of  the 
transition  matrix,  qAC: 

E[uf]  =  ((-l)2  (2i>+c)+(l)2  (2b+c))  (26) 

=  (4fc+2c) 

Since  the  Uj  are  independent  and  identically  distributed  random  variables,  by  the 
Central  Limit  Theorem,  Xn  is  approximately  normally  distributed  with  mean 
/ix=E[Xn]=0  and  variance  a2x=Var[Xn]=n(Var[Uj])=n(4b+2c),  for  large  n 
[Ref.  5:  p.  289].  The  same  derivation  can  be  applied  to  the  vertical  position 


r.  *  £ 


(27) 


where  v,  is  the  vertical  movement  at  time  i. 

Therefore,  (Xn,Yn)  is  approximately  bivariate  normal  with  parameters 
/z=jzx=jxv=0,  and  a2=a2x=a2y=n(4b+2c).  Furthermore,  the  covariance: 


Cov[Xn,yj  =  e\xh  y„  ]  -  E[x„  ]E[yh  ]  (28) 

-Wmrm\ 

where  the  (XnYn)e(-n2,...,-l,0,l,...,n2)  with  probability  distribution 
[p(-n2),...,p(-l),p(0),p(l),...,p(n2)].  Because  of  the  symmetry  of  qAC,  we  can  see  that: 
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(29) 


pi-n2)  =  PiX^-nJ^nhPiX^nJ^-n) 
=  P(Xn=n,YH=n)+P(XH=  -n,YH=  -rt) 

=  Pin2) 


Thus: 

E[XJJ  -  in2)pin2M-n2)p(-n2)+...+ 

(l>p(l)+(-  1)/K- D+(0)/K0)  (30) 

=  0 

Therefore,  Cov[Xn,Yn]=0  so  Xn  and  Yn  are  uncorrelated  [Ref.  5:  p.310].  The  result 
is  an  approximately  circular  normal  distribution  with  mean  0  and  variance  n[4b+2c] 
when  n  is  large.  Note  that  this  is  true  independent  of  the  values  of  a,  b  and  c.  The 
important  fact  is  that  the  transition  matrix  maintains  horizontal  and  vertical 
symmetry.  For  most  central  limit  theorem  applications,  large  n  is  taken  to  be  greater 
than  20.  As  an  example,  for  a  target  speed  of  15  knots,  a  mean  time  on  leg  of  one 
hour  and  a  cell  size  of  5  nautical  miles,  we  find  from  Equations  (13)  and  (14)  that 
just  over  two  hours  would  be  required  to  meet  this  central  limit  theorem  requirement. 

We  can  now  relate  the  number  of  iterations  to  the  radius  containing  86 
percent  of  the  target  density,  R86.  The  radius  R=(X2n+Y2n)'/2.  Assuming  the 
requirements  of  the  central  limit  theorem  are  met,  X„  and  Yn  are  Normally 
distributed,  and  because  the  are  uncorrelated,  they  are  also  independent.  Thus  XJo 
and  Yn/a  are  independent  standard  normal  random  variables,  and: 
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(31) 


x2  1?  Xl*  1? 


is  approximately  a  Chi-Squared  random  variable  with  two  degrees  of  freedom.  As 
discussed  earlier,  86  percent  of  the  density  is  contained  within  two  standard 
deviations,  we  will  call  this  with  the  subscript  TM  representing  the  Transition 

Matrix  method.  Let  W=(cell  width),  then  for  large  n: 

R862tm  =  (2  a)2 

=  (32) 

=  n(16b+4c)W2 


where  ^£6^  is  measured  in  nautical  miles.  For  the  values  of  a,  b  and  c  used  in  qAC, 
this  results  in  the  relationship: 

R862  =  2.54  nW2  (33) 

This  is  very  close  to  the  experimentally  obtained  relationship  of  Equation  (14)  used 
in  ASWTDA. 

In  ASWTDA,  the  Transition  Matrix  method  is  matched  to  the  random  tour 
by  equating  Equation  (32)  with  2E[R2T]  from  Equation  (16),  even  when  n  is  small. 
The  method  keeps  a,  b  and  c  constant  while  making  n  be  some  integer  other  than  the 
number  of  information  updates  (recall  that  these  updates  occur  at  regular  intervals 
of  A  =6  Min.).  It  would  be  more  convenient  to  make  n  the  same  as  the  number  of 
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information  updates.  If  V  is  not  too  large  compared  to  cell  size,  this  can  be  achieved 
using  the  exact  expression  of  Equation  (24)  instead  of  the  approximate  expression  of 
Equation  (13)  by  letting  a,  b  and  c  depend  on  time.  The  crucial  question  is  whether 
the  increment  (16b+4c)W2  from  Equation  (32)  can  be  made  as  large  as  the  increase 
of  Equation  (24)  in  a  period  of  length  A.  The  slope  of  Equation  (24)  is  at  most  4V2X 
and  the  maximum  value  of  (16b+4c)  is  41,  so  the  crucial  inequality  is  4V2XA<4W2. 
Letting  A =0.1  hour,  we  obtain  the  requirement: 

V  <;  W.  —  (34) 

N 

where  V  is  measured  in  knots. 

Note  that  for  T>2A,  the  slope  of  Equation  (24)  is  a  constant  value  of 
4V2A,  therefore,  a,  b  and  c  may  be  constant  for  T>2A. 

One  can  see  from  Equation  (34)  that  there  would  be  only  few  instances 
where  a  greater  number  of  diffusions  would  be  necessary.  For  example,  for  a  mean 
time  on  leg  of  one  hour  and  a  cell  size  of  five  nautical  miles,  V<15.8  Kts.  These 
instances  could  be  handled  on  a  case  by  case  basis,  and  n  could  still  be  made  a 
convenient  integer  value. 


1  This  occurs  when  b=0.25,  i.e.,  all  probability  in  the  comer  cells.  Though  it  is 
recognized  that  this  would  result  in  some  cells  being  inaccessible,  it  is  the  limiting 


2.  Primary  Target  Diffusion 

In  ASWTDA,  the  transition  matrices  are  such  that  the  target  may  move 
only  as  far  as  one  cell  with  each  application  of  q.  This  does  not  result  in  errors  for 
the  Area  Clearance  target,  but  the  method  of  calculating  the  time  interval,  dt,  for 
Primary  Target  diffusion,  may  result  in  some  inaccuracy.  For  example,  the  Primary 
Target  may  be  going  faster  than  the  average  speed,  and  therefore  may  travel  farther 
than  one  cell  length  in  a  time  step  dt,  as  illustrated  in  time  step  number  8  in  Figure  4. 
The  transition  for  this  time  step  can  either  be  ignored,  or,  as  is  done  in  ASWTDA, 
a  count  can  be  placed  in  the  closest  adjacent  cell,  inaccurate  in  either  case.  Thus 
there  is  a  tradeoff  involved  in  the  choice  of  dt.  Too  large  a  dt  results  in  the  target 
being  able  to  go  farther  than  one  cell  length,  resulting  in  the  situation  just  described; 
too  small  a  dt  results  in  too  few  transitions  to  adjacent  cells  resulting  in  a  high 
probability  of  remaining  in  the  center  cell  of  the  transition  matrix  thus  requiring  more 
frequent  application  of  the  transition  matrix,  increasing  computation  time.  In  order 
to  prevent  a  target  from  traveling  farther  than  one  cell  length,  the  appropriate  time 
step  would  be: 

dt  =  CeU  Length  (35) 

Maximum  Target  Speed 

It  is  not  intended  here  to  state  an  optimum  method  of  choosing  dt,  simply 
that  diffusion  accuracy  is  dependent  on  this  choice,  which  in  turn  depends  on  the 
target  velocity  distribution,  and,  as  shown  by  the  need  for  the  cosine  correction  term 
in  certain  cases,  on  average  target  course  as  well.  Note  however  that  the  average 


31 


target  course  has  no  meaning  when  the  variance  in  target  course  is  large,  with  the 
extreme  case  being  a  randomly  touring  target.  This  aspect  of  the  problem  warrants 
further  study.  Of  course,  another  way  to  prevent  this  error  is  to  expand  the  transition 
matrices,  q,  to  say  a  5x5  matrix.  The  tradeoff  here  would  be  an  increase  in  both 
computation  time  and  storage  requirements. 
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VI.  TARGET  COMBINATION  AND  CREDIBILITY 

A.  TARGET  COMBINATION 

The  question  arises  as  to  what  specifically  is  described  by  the  Area  Clearance 
and  Primary  Target  density  matrices  and  the  operations  performed.  Two  elementary 
models  are  presented  here,  these  are:  1.  A  single  target  which  may  be  of  one  of  two 
types,  2.  Two  separate  targets.  An  attempt  will  be  made  in  Section  3  to  match  each 
of  these  models  with  the  operations  performed  in  ASWTDA.  This  is  not  intended 
to  represent  all  possible  target  models,  merely  two  that  were  chosen,  that  have  a 
sound  probabilistic  basis,  to  attempt  to  understand  the  operations  of  ASWTDA. 

1.  Single  Target 

In  rhis  scenario,  it  is  known  that  there  is  a  single  target  within  the  analysis 
area  at  time  t=0,  but  it  is  not  known  of  what  type.  Assume  the  target  is  of  one  of 
two  types,  a  or  b,  and  we  presumably  have  some  knowledge  of  the  a  priori  probability 
distribution  of  target  type.  Furthermore,  it  is  known  that  the  type  a  target  remains 
within  the  analysis  area,  and  the  type  b  target  may  leave  the  area  and  not  reenter. 
For  this  scenario: 
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(36) 


Au(i)  =  PiEJt)  and  a|ff,) 

A'/t)  =  PiE^it)  and  a  \  Htl(t)) 

B,/t)  =  P(£o<0  and  b  |  Ht) 

B,>)  =  PO^tt  and  b|iy,/(0) 

C(r)  =  P((outside  S)  and  b  |  Ht) 

C*(t )  =  P((outside  S)  and  bj  HtI(t)) 

so  |  A(t)  |  =P(a  |  Ht),  likewise  for  target  type  b,  and  |  A(t)+B(t)  |  +C(t)=  1. 

There  are  two  separate  motion  models  involved  described  by  the  following 
transition  matrices: 


Qa(i,j\m,ri)  =  P(Eid(t+ 1)  and  a\EnA(t)  and  a) 
Q\i,j\m,n)  =  P(EtJ(t+l)  and  b\Em/t)  and  b) 


Thus: 


^(r+l)  =A\t)*Q°  (38) 

B(r+ 1)  =  B '(*)*<?* 

and  C(t+ 1)= 1- 1  A(t+ 1)+B(t+ 1)  | . 

These  matrices  may  be  updated  for  received  information  using  Bayes’ 
theorem  as  follows: 


P(£(J(r)a  | //,/(*)) 


Pm\EiJ(t)aHl)P(Eij(t)a\Ht) 

pm\Ht ) 


P(Eifi)  b  1  Htl(t)) 
P((outside  S)  b  |  HtI(t )) 


P(/(Q  |  E^bH)  PjE^t)  b  |  Ht)  (39) 

pmw 

P(/(f)  |  (outside  S)b//,)  P((outside  5)  b  j  H) 
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As  discussed  in  Section  II.C.l,  we  would  expect  to  assume  that 
P(I(t)  j  Ejj(t)Hta)  and  P(I(t)  |  Ejj(t)H,b)  are  both  independent  of  H„  but  the 


nondetection  probability  may  still,  in  principle,  depend  on  target  type. 
To  simplify  notation,  the  following  are  defined: 


=  Pm\EJt)&Ht) 
btp)  =  PW)  \Ei/t)bHl) 
c(t)  =  P(l(t)  |  (outside  S)bHt) 

Thus,  Equations  (39)  become: 


= 

b;p)  - 

CTO  = 


pm  w 
Kmp) 
pm  \h,) 

c(QC(Q 

p(m\Ht) 


(40) 


(41) 


and  since: 


pm \h,)  =  E  [flyco^o  +  tjL  +  c(oc(o  (42) 

iJeS 

we  can  see  that  |  A+(t)+B+(t)  |  +C+(t)=l. 

2.  Two  Separate  Targets 

In  this  model,  there  are  two  targets,  one,  a  type  a  target  is  in  the  analysis 
area,  and  the  other,  type  b,  possibly  outside  the  analysis  area.  For  this  model,  the 
following  are  defined: 


AJt)  =  P(E./t)  la/?,) 

A>)  -  PiE^lzH'Kt)) 

BJ&  *  PVJ&W  (43) 

2&0  =  P(EJt)\bHtm 
C(t)  =  P(outside  S  |b/f,) 

C'(t)  =  P(outsideS|bff,/(*)) 

Motion  updating  is  equivalent  to  that  for  the  single  target  model. 

The  density  matrices  are  updated  for  received  information  using  Bayes’ 
theorem  as  follows: 


BlU) 

C\t) 


Pjlit)  lE^&HJA'p) 

Y,Pm\EJf>*Ht)AJt> 

iJeS 

_ Pm\EtJ(t)bHt)BJt) _ 

£  [PtMlEdtobHJBJto]  *  Pm  I  (outside  S)bHt)  C(i) 

ijeS 


(44) 


3.  Analysis 

a.  Single  Target  Model 

To  compare  the  operations  performed  in  ASWTDA  with  the  single 
target  model  as  described  above,  we  must  make  an  assumption  about  the  meaning 
of  the  target  density  matrices.  Let  Tac  and  Tpt  be  the  events  the  target  is  of  type 
Area  Clearance  and  Primary  Target,  respectively.  We  will  assume  the  following: 


AC,p  =  PiEJQT^H,)  (45) 

PT,/»  -  nEJi)T„\Ht) 

Note  however  from  Equations  (y)  and  (11)  that  following  an  information  update, 
|  AC+(t)  |  =  j  PT+(t)  |  =1.  ASWTDA  then  performs  an  additional  operation  to  form 
what  we  will  call  AC++(t)  and  PT++(t),  which  is  a  renormalization  of  the  density 
matrices  such  that: 


AC+'{t)  =  AC*(t)\AC(t)\ 
FT*(t)  =  FT (t)  \  PT(t)  | 

In  the  notation  of  the  previous  section,  this  is  equivalent  to: 


a.  (t)A  (t) 


(46) 


a(t)A(t) 

+  b.  (t)B.  (t) 

B  :(t)  =  -'ll..— f  A :  |  B(l) 
|  b(t)B(t)  I  1  w 


(47) 


These  equations  must  be  compared  with  the  results  of  Equations  (41). 
Using  the  very  simplified,  two  cell  example  shown  in  Figure  5,  these  are  clearly  not 
equivalent  calculations.  Furthermore,  although  Reference  1  states  that  the  Primary 
Target  may  leave  the  analysis  area,  following  the  target  motion  update,  the  Area 
Clearance  matrix  is  again  renormalized  such  that  |  AC  |  =  1- 1 PT  | .  This  is  a  deviation 
from  the  single  target  model  in  that,  following  each  motion  update,  C(t)=0.  So, 
rather  than  leave  the  area,  the  target  type  shifts  to  an  Area  Clearance  target. 

ASWTDA  also  handles  positive  information  differently  than  negative 
information.  As  discussed  in  Chapter  IV,  PT(t)=0  for  t<t0.  Thus,  the  Primary 
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Target  distribution  is  based  on 
Area  Clearance  type  motion  for 
t<T0,  and  Primary  target  motion 
for  t>t0.  This  too  is  a  deviation 
from  the  single  target  model  and 
results  in  an  inaccurate  Primary 
Target  distribution  to  a  degree 
that  is  dependent  upon  the 
difference  in  the  motion  models 
and  the  elapsed  time  from  t=0 
to  t=t0. 

As  can  be  seen 

from  Equations  (39),  both  positive  and  negative  information  can  be  handled  equally 
well  using  Bayes’  theorem,  provided  the  appropriate  probabilities  are  known. 

It  is  conceivable  that  receipt  of  positive  information  may  alter  the 
distribution  of  target  type.  This  is  properly  handled  by  Equations  (39).  For  example, 
if  we  make  the  assumption  that  some  received  positive  information,  I(t),  could  not 
possibly  been  caused  by  target  type  a,  in  Equations  (39),  P(I(t)  |  Ejj(t)aH,)=0  for  all 
(i,j),  thus  P(Ejj(t)  |  H,I(t))=0  for  all  (i,j),  and  all  probability  is  assigned  to  target  type 
b.  In  order  to  perform  this  operation  properly  in  ASWTDA,  it  is  necessary  that  the 
density  matrix  for  the  Primary  Target  be  nonzero  for  t<to,  otherwise  P(I(t)  |  H,)=0 
which  results  in  division  by  zero. 


Figure  S  Comparison  of  Bayesian  Update  and 


ASWTDA  Update  Methods 


38 


b.  Two  Target  Model 


To  compare  ASWTDA  with  the  two  target  model,  we  must  assume 
the  target  density  matrices  are  defined  as  follows: 

ACJf)  =  P(EJ')  |  TACHt) 

AC‘/I)  =  P(£,/0 1  TicH,m  4g) 

PTJt)  =  P(EJt)\T„H,) 

PT’Jt)  =  P(EJt)  |  T„H,m 

The  two  target  model  clearly  is  not  what  is  intended  in  ASWTDA, 
because,  as  stated  in  Reference  1,  "...each  of  the  density  displays  represents  position 
information  for  the  same  single  target  [Ref.  1:  p.  9]".  This  is  further  reinforced  by 
the  fact  that,  following  each  information  and  motion  update,  ASWTDA  adjusts  the 
density  matrices  such  that  |  AC  |  +  |  PT  |  =  1. 

4.  Modifications  to  Comply  With  a  One  or  Two  Target  Model 

Changes  may  be  made  to  ASWTDA  so  that  the  operations  conform  to 
either  of  the  two  motion  models  described  above. 

In  order  to  conform  to  the  single  target  model,  the  following  changes  must 

be  made: 


1.  Maintain  a  Primary  Target  density  matrix  from  time  t=0,  which  follows  the 
motion  model  that  is  believed  to  describe  the  Primary  Target.  If  the  motion 
model  for  the  Primary  Target  is  unknown,  multiple  possible  scenarios  could 
be  initialized  at  t=0.  Then,  at  t=to,  the  probability  of  each  of  these  scenarios 
being  the  true  motion  model  can  be  assessed.  Another  option  is  that  upon 
receipt  of  positive  information,  a  Primary  Target  motion  model  may  be 
developed  and  started  from  t=0  to  calculate  a  density  PT(t0). 


39 


2.  Initialize  the  Area  Clearance  and  Primary  Target  prior  distributions 
independently  to  what  is  believed,  a  priori,  about  the  individual  targets.  If  it 
is  known  that  there  is  one  target  in  the  analysis  area,  then 
|  AC(0)  |  +  |  PT(0)  |  =  1. 

3.  If  the  definitions  given  in  the  single  target  analysis  section  above  are  what  is 
intended  in  ASWTDA,  then  updating  for  information  should  be  performed 
using  Equations  (41). 

4.  As  the  Primary  Target  diffuses  out  of  the  analysis  area,  the  Area  Clearance 
matrix  should  not  be  renormalized  to  include  this  lost  target  density. 

5.  Since  the  operator  is  concerned  about  a  single  target,  his  attention  should  be 
focused  on  the  combined  target  density,  which  is  the  best  knowledge  available 
on  the  single  target. 


If  the  intention  in  ASWTDA  is  that  there  are  two  targets,  then  in  order  to 
conform  to  the  two  target  model,  in  addition  to  number  1  above,  the  following 
modifications  must  be  made: 


1.  Initialize  the  Area  Clearance  and  Primary  Target  prior  distributions 
independently  to  what  is  believed,  a  priori,  about  the  individual  targets. 
Following  the  two  target  model,  |  AC(0)  j  =1,  and  |  PT(0)  |  <1. 

2.  As  the  Primary  Target  diffuses  out  of  the  analysis  area,  only  the  Primary 
Target  density  matrix  should  be  affected.  Therefore  |PT(t+l)|<|PT+(t)| 
whereas  |  AC(t+ 1)  |  =  |  AC+(t)  |  =  1. 

3.  Update  the  target  density  matrices  for  information  using  Bayes’  theorem  of 
the  form  of  Equations  (44). 

4.  The  density  matrices  should  not  be  combined  into  a  single  target  density 
matrix  as  they  represent  two  separate  targets. 
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B.  CREDIBILITY 


1.  A  Proper  Application  of  Credibility 

A  decision  aid  must  account  for  the  fact  that  the  emissions  which  resulted 
in  the  positive  information  report  may  have  come  from  a  target  different  form  the 
target  of  interest.  This  may  be  as  simple  as  assuming  that  it  does  not  occur.  Let  T 
be  the  event  the  emissions  are  caused  by  the  target  of  interest,  let  I(t)  be  the  event 
a  positive  report  is  made,  and  assume  the  following  data  are  known: 

ccu(t)  =  P(I(0T\EJt))  for  all  (iJ)eS 
P JO  =  PmTc\EJt))  for  all  ( iJ)eS  (49) 

Y(r)  =  P(/(r)  T  |  outside  S) 

T](r)  =  P(I(t)Tc  | outsider) 

Thus  is  the  event  that  a  report  is  made  and  the  target  of  interest  is  not  the 
cause,  and  a1J(t)+)Siij(t)=P(I(t)  |  Eitj(t)),  where  P(I(t)  |  Ejj(t))  is  discussed  earlier  in 
Section  II.C.l.  By  the  definition  of  conditional  probability: 


P(E./t)  |  /(/))= 


PmTlE'/tVPjEJt)) 

pm)) 


+  P(7(f)rclE./r))P(E,>>(r)) 

P(/(0) 


(50) 


Note  that  P(I(t))  can  be  written  as  follows: 


pm))  = 


pmm 

P(T\M) 

pm)Tc) 

p(Tc\m 


(51) 


Letting  pij(t)=P(Ejj(t))  and  P(t)  be  the  matrix  with  elements  pi(j(t),  and  letting  a(t) 
be  the  matrix  with  elements  ceij(t)  and  /3(t)  the  matrix  with  elements  ^(t),  Equation 


(50)  can  be  written: 
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rout  MO  ■  15MBM  ,(r  |j(i» 

“  «M  mtoro  ' 

_ _ W|/<0)-( _ MML- 


|a(0mhY«)(l-|W|) 


PWmi^Wd-IWI) 


h 


(52) 
1 7(0) 


Let  P(T  |  I(t))  be  an  input  called  credibility,  C.  If  we  assume  >9id(t)=>3  for 
all  (i,j)  and  all  t,  i.e.,  if  I(t)  is  not  caused  by  the  target  of  interest,  it  does  not  matter 
where  it  is  located,  and  if  we  make  the  further  assumption  that  y(t)=T7(t)=0  for  all 
(t),  i.e.,  we  can  not  possibly  receive  a  contact  report  if  the  target  is  outside  S,  even 
if  the  target  is  not  the  cause  of  the  report,  Equation  (52)  becomes  [Ref.  6:  pp.  8-9]: 


P(Ei/t)\I(t))  = 


C  +  p.p)(\-C) 


(53) 


In  other  words,  the  posterior  probabilities  are  an  average  of  the  prior  probabilities 
(with  probability  (1-C)),  and  the  posterior  probabilities  that  would  hold  when  C=1 
(with  probability  C). 


2.  Credibility  in  ASWTDA 

It  is  difficult  to  determine  what  is  implied  by  the  use  of  the  confidence 
factor,  c/,  in  ASWTDA.  This  is  in  part  due  to  the  fact  that  it  is  not  clearly  defined 
what  is  described  by  the  Area  Clearance  and  Primary  Target  density  matrices.  An 
attempt  will  be  made  to  fit  the  above  model,  however,  based  on  the  single  target 
model. 

Let  T  be  the  event  the  emissions  which  resulted  in  the  positive  report  were 
caused  by  the  Primary  Target,  and  I(t0)  the  event  that  positive  information  is  received 
at  time  (t0).  One  infers  from  Reference  1  that  the  confidence  factor,  cf=?(T\  I(t0)), 
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i.e.,  cf  is  the  probability  that  the  received  emissions  were  caused  by  the  Primary 
Target  given  the  receipt  of  the  positive  information  at  time  (t0).  For  the  analysis,  the 
following  assumptions  will  be  made: 


1-  -ffZ,ij=p(Kto)T  I  Ei(j)  for  all  (ij)eS. 

2.  C=P(T  |  I^)),  i.e.,  C  is  the  probability  the  emissions  which  resulted  in  the 
positive  information  report  were  caused  by  the  Primary  Target,  given  the 
receipt  of  positive  information  at  time  (to). 

3.  P(I(t0)Tc  |  Ejj(t0))is  constant  for  all  (i,j). 

4.  P(I(t0)T  j  outside  5)=P(I(t0)Tc  |  outside  S)=0 


Equation  (53)  then  becomes: 


/,(£ls/{r0)  (Primary  Target)  \IitJ)  =  P'/tJC  +  ^(^(1  -  C)  (54> 
This  is  equivalent  to  Equation  (10)  with  two  exceptions: 

1.  The  multiplication  term,  pr,  used  in  ASWTDA  is  the  confidence  factor 
reduced  further  by  excluding  that  portion  of  B°  which  falls  outside  the  analysis 
area  or  on  land.  Thus  the  confidence  of  the  report  is  dependent  upon  the 
location  of  the  analysis  area  in  relation  to  the  positive  report. 

2.  The  second  term  of  Equation  (54)  is  omitted  altogether. 

Thus  we  can  see  that,  even  given  these  assumptions  about  the  target  and 
the  received  information,  the  use  of  the  confidence  in  ASWTDA  falls  short  of  a 
meaningful  application  of  target  credibility.  This  may  be  easily  fixed  by  utilizing 
Equation  (54)  in  ASWTDA,  use  of  which  would  correspond  to  making  assumptions 
1,  2,  3  and  4  above. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  summary,  there  are  four  areas  which  need  to  be  addressed  in  the  continuing 
development  of  ASWTDA  to  prevent  possibly  erroneous  results.  These  are: 

1.  The  number  of  targets  in  the  system 

2.  Proper  incorporation  of  positive  information 

3.  Primary  target  motion  prior  to  receipt  of  positive  information 

4.  The  proper  use  of  information  credibility 

The  question  of  how  many  targets  there  are  and  what  is  meant  by  an  Area 
Clearance  target  or  a  Primary  Target  has  a  profound  affect  on  how  the  other  points 
are  handled.  This  issue  must  be  resolved  prior  to  addressing  further  issues.  Once 
it  is  decided  how  many  targets  are  being  manipulated,  it  can  be  decided  how  the 
target  or  targets  are  to  move,  and  to  which  target  or  targets  positive  information  is 
applied.  Also,  it  must  be  decided  what  is  meant  by  the  confidence  factor  in  relation 
to  positive  information,  and  to  apply  it  properly. 

Other  aspects  which  may  not  significantly  degrade  the  results  obtained  from 
ASWTDA,  but  may  affect  the  speed  or  simplicity  of  computation  and  overall  solution 
accuracy  are  use  of  the  exact  expression  for  R86  for  a  randomly  touring  target  and 
the  adjusting  of  the  values  of  a,  b  and  c  in  the  Area  Clearance  transition  matrix  so 
that  the  diffusion  of  the  Area  Clearance  target  occurs  at  regular  convenient  intervals. 
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